Preprint typeset in JHEP style - HYPER VERSION 



pre-print xxx 



Determination of 77 ^ tt+tt tt^ Dalitz plot slopes and 
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Abstract: We have studied, with the KLOE detector at the DA$NE <I>-Factory, the 
dynamics of the decay tj 7r~^TT~Tr^ using rj mesons from the decay — > 777 for an 
integrated luminosity C = 450 pb^^. Prom a fit to the DaUtz plot density distribution 
we obtain a precise measurement of the slope parameters. An alternative parametrization 
relates the tt^tt^tt^ slopes to that for rj 3tt^ showing the consistency of KLOE results 
for both channels. We also obtain the best confirmation of the C-invariance in the rj — > 
7r+7r~7r° decay. 
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1. Introduction 

The decay rj ^ 3tt violates iso-spin invariance. Electromagnetic contributions to the pro- 
cess are very small IH and the decay is induced dominantly by the strong interaction via 
the u, d mass difference. The 77 ^ Svr decay is therefore an ideal laboratory for testing 
chiral perturbation theory, ChPT. A three body decay^ is fully described by two variables. 
We can choose two of the pion energies {E^, Eq) in the rj rest frame, two of the three 

two pion masses squared (m^ , m^Q; ^^0+) ^^^^ called {s,t,u). Note that E^ is linear in 

m^Q and so on, cyclically. We use the Dalitz variables, X, Y which are linear combinations 
of the pion energies: 



a.i) 



where Q is the decay "Q-value". The decay amplitude is given in Q as: 

where = [ml — fh'^)/{m^ — m^) and m = {niu + md)/2 is the average u, d quark mass. 
-Ftt = 92.4 MeV is the pion decay constant and 27l(s, t, u) must come from theory. From eq. 



""^Both T) and tt are spinless, therefore there is no preferred direction. 
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1.2 it follows that the decay rate for rj ■k'^tt -k^ is proportional to A ^. The transition 
r/ — > Svr is therefore very sensitive to A if the amplitude 3Jt were known. At lowest order 
in ChPT: 

m{s,t,u) = —^ — ^. (1.3) 



From eq. |l.3| , 0] one finds {q tt'^tt vr ) =66 eV to be compared with the measured 



width of 295 it 16 eV y]. A one-loop calculation within conventional chiral perturbation 
theory (ChPT) improves considerably the prediction: 



r""^" (r/ ^ 7r+7r"7r° ) ~ 167 ± 50 eV. (1.4) 

but is still far from the experimental value. Higher order corrections ||5| help but do not 
yet bring agreement with measurements of both total rate and Dalitz plot slopes. Good 
agreement is found combining ChPT with a non perturbative coupled channels approach 
using the Bethe Salpeter equation Q. 

Therefore a precision study of the r] ^ 3tt Dalitz plot, DP, is highly desirable. The 
amplitude squared is expanded around X = y = in power of X and Y 

\A{X,Y)\'^ 1 + aV + bV"^ + cX + dX^ + eXY + .... (1.5) 

The parameters (a, b, c, d, e, ... ) can be obtained from a fit to the observed DP density and 
should be computed by the theory. Any odd power of X in A{X, Y) implies violation of 
charge conjugation. 



2. The KLOE detector 

Data were collected with the KLOE detector at DA<I>NE the Frascati e^e~ collider, 
which operates at a center of mass energy W = ~ 1020 MeV . The electron and positron 
beams collide with a crossing angle of vr — 25 mrad, resulting in a small momentum [p^ ~ 
13 MeV/c in the horizontal plane) of the produced (j) mesons. The KLOE detector consists 
of a large cylindrical drift chamber (DC), surrounded by a fine sampling lead-scintillating 
fibers electromagnetic calorimeter (EMC) inserted in a 0.52 T magnetic field. 

The DC 4 m diameter and 3.3 m long, has full stereo geometry and operates with a 
gas mixture of 90% helium and 10% isobutane. Momentum resolution is cr{p±)/p± < 0.4%. 
Position resolution in r — is 150 ^m and ~ 2 mm. Charged tracks vertices are 
reconstructed with an accuracy of ~ 3 mm. 

The EMC is divided into a barrel and two endcaps, for a total of 88 modules, 
and covers 98% of the solid angle. Arrival times of particles and space positions of the 
energy deposits are obtained from the signals collected at the two ends of the calorimeter 
modules, with a granularity of ~(4.4 x 4.4) cm^, for a total of 2240 cells arranged in five 
layers. Cells close in time and space are grouped into a calorimeter cluster. The cluster 
energy E is the sum of the cell energies, while the cluster time t and its position r are 
energy weighted averages. The respective resolutions are cte/E = 5.7%/ E (GeV) and 
at = 57 ps/^/E (GeV) 100 ps. 
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The KLOE trigger |10| is based on the coincidence of at least two energy deposits in 
the EMC above a threshold that ranges between 50 and 150 MeV. In order to reduce the 
trigger rate due to cosmic rays crossing the detector, events with a large energy release in 
the outermost calorimeter planes are vetoed. 



3. Signal selection and efficiency 

This analysis refers to ~ 450 pb^^ collected at DA$NE in years 2001/02 corresponding to 
~ 1.4 10^ (j) mesons produced. 

At KLOE 7] mesons are produced through the radiative decay — > 777 . Accounting 
for the product of BR's: BR{(p rj'j ) xBR(t/ 7r+7r~7r° ) ~ 2.9 x 10^'^ we expect about 
four millions of ry — > tt^tt^tt^ events. A larger Monte Carlo (MC) sample corresponding to 
about 5 times the amount of data has been used to study efficiencies and backgrounds. 

Note that the recoil photon is almost monochromatic, with rec ~ 363 MeV, well 
separated from the softer photons from vr'^ decay. 

A photon is defined as an EMC cluster not associated to a DC track. We further 
require that \{t — r/c)\ < 5at, where t is the arrival time at the EMC, r is the distance of 
the cluster from interaction point, IP, c is speed of light. The events selection is performed 
through the following steps: 

1. Events are first selected by a very loose offline filter to remove machine background 
(FILFO) and an event selection procedure (EVCL) assigning events into categories 
0- 

2. We then require two opposite curvature tracks intersecting at a point (vertex) inside 
a cylinder with r < 4 cm, jz] < 8 cm centered at the IP. We require also three photons 
with 21° < 9^ < 159° and > 10 MeV. The angle between any photon pair must 
be > 18° to remove split showers. 

3. < 800 MeV. 

4. A constrained kinematic fit is performed imposing 4-momentum conservation and t = 
r/c for each photon. We retain events with a probability P{x^) > 1%; corresponding 
to < 18. The fit significantly improves the photon energy resolution. The 
distribution is in reasonable agreement with MC prediction, as shown in fig.|l| varying 
the cut on P{x^) in the range [ 0.01, 0.15] (x^ < 18tox^ < 10) has no significant 
effect on the analysis results, see section |5[ 

5. Finally we require: 

(a) 320 MeV< E.y^ rec < 400 MeV for the recoil photon, to reduce residual back- 
ground from (j) Ks Kl events. 

(b) E^+ +E^~ <550 MeV, to reduce residual background from (p — > vr^vr^vr'^ events. 

(c) ?7i(77) for the two softest photons must satisfy 110< m^^ <160 MeV, to reduce 
residual background from r/ — > 7r"'"7r~7r'' decays with vr'^ — > e~^e~'y and from 
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Figure 1: distribution for the kinematic fit. Left: linear scale. Right: log scale. 



• LOTT^ with Lo — > vr^vr vr'' ; and to eliminate the residual background from 
777 events with rj — > tt^tt^^. 



The selection efficiency is determined with the MC program [11 1 and checked with data 
control samples. In particular: 



1. The trigger efficiency evaluated by MC is 99. 
for the trigger sectors multiplicities. 



with excellent data-MC agreement 



2. The effects of EVCL and FILFO are evaluated using a downscaled set of non filtered 
data with less stringent cuts in order to get a "minimum bias" sample. On signal 
events the efficiency of the minimum bias selection is 99.88%. We have found that the 
EVCL procedure introduces a signal loss of ~1.5%, as expected also from MC. The 
corresponding bias on the Dalitz plot parameters has been included in the systematic 
error. No bias is introduced by the FILFO procedure. 

3. The tracking and vertexing efficiencies have been estimated from the data-MC ratio 
observed in a sample of ^ ir^-K'^Tr^ events with charged pion momenta in the same 
range as those from the r] — > tt^tt^-k^ decay [12|. These events can be selected with 
low background requiring the detection of the photons associated to the vr'^ in the 
EMC and only one track in the DC and thus are suited to study on data the single 
charged track reconstruction efficiency, etrki and the charged vertex reconstruction 
efficiency, evtx: e?rk^vtx, data/etrk^vtx, MC = 0.974 ± 0.006. This ratio is constant for 
all momenta, introducing no bias in the Dalitz plot distribution. All variables used in 
the fit are evaluated in the 77 rest frame, which in the laboratory has a momentum of 
~363 MeV. Therefore to each momentum bin in the rest frame corresponds a wider 
interval in the lab; MC-data discrepancies are further diluted by this effect. 

4. A correction to the MC detection efficiency for low energy photons has been obtained 
by comparing the photon energy spectrum of a data subsample to the expected MC 
spectrum; the average correction factor is 0.964. 

The overall selection efficiency, taking into account all the data-MC corrections is found 
to be e = (33.4 it 0.2)%. The expected background contamination, obtained from MC 
simulation is 0.3%. 
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After background subtraction we remain with 1.34 • 10 events. 
The DaUtz plot density is shown in fig.||. 




Figure 2: DP distribution for the whole data sample. The plot contains 1.34 millions of events in 
256 bins. 

The signal selection efficiency e {X, Y) as function of the DP point is obtained by MC, 
for each {X, Y) bin, as the ratio: 



e{X,Y) 



Nrec iX,Y) 
iVgen {X,Y) 



(3.1) 



where A'rec, gcn(-'^; ^) are respectively the reconstructed and generated DP populations. 
This approach accounts for resolution effects as long as the MC correctly reproduces the 
Dalitz plot shape; a first estimate of the Dalitz plot parameters to be used in the final 
MC was obtained from a preliminary fit to a data subsample. The efficiency e {X, Y) has 
a smooth behavior all over the entire DP. The projections of e{X, Y) are shown in figjsj. 
While the efficiency appears to be rather flat on X (and symmetric as expected), it decreases 
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Figure 3: Left: Efficiency vs X. Right: Efhciency vs Y. 

approximately linearly with Y. In fact a large Y value means a low-momentum vr^ in the 
decay to which corresponds a lower tracking/vertexing efficiency. The resolutions from MC 
on the DP variables {X, Y) are shown in fig.^. The Y variable, which is proportional to 
the vr*^ kinetic energy, is evaluated |13| as the average between the "direct" determination 
obtained from the energy and direction of the two clusters associated to the vr" 77 
decay and the "indirect" determination : Tq = — {E^+ + E^- ) — M^o . Due to our 
excellent momentum resolution for charged tracks, the core of both distributions in figj^ 
can be fitted with a gaussian with a = 0.02. 
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Figure 4: X (left) and Y (right) resolution from MC. 

4. Fit of Dalitz plot 

The expected Dahtz density is taken as: 

r(X, Y) = \A{X, = N{1 + aY + bV^ + cX + dX^ + eXY + ...). (4.1) 

with N being a normahzation constant. The fit to the Dahtz plot is done in two dimensions, 
minimizing the function. Bins intersecting the Dalitz plot boundary are not included 
in the fit. The fit procedure has been tested on MC by verifying that the fit reproduces in 
output the same input values of the DP parameters. 

The fit results for different forms of \A\'^ and for AX=Ay=0.125 (154 bins fitted) are 
shown in Table ||. A fit with only quadratic terms gives a very low C.L. of 0(10"^) or less. 
Including cubic terms as 

T{X,Y) = N{l+aY+bY'^+cX+dX'^ + eXY+fY^+gX^ + hX'^Y+lXY'^) (4.2) 

results in much better fits with C.L. > 70 %. In particular the coefficients / of the Y^ 
term and d of the term are clearly required while the other ones {g, h, I) turn out to be 
consistent with zero. As expected from C-invariance c and e are consistent with zero. 



dof 


CL 


axlO^ 


6x10^ 


cxlO^ 


dxlO^ 


exlO^ 


/xlO^ 


147 


73% 


-1090±5 


124±6 


2±3 


57±6 


-6±7 


140±10 


149 


74% 


-1090±5 


124±6 




57±6 




140±10 


150 


< 10-6 


-1069±5 


104±5 








130±10 


150 


< 10-8 


-1041±3 


145±6 




50±6 






151 


< 10-6 


-1026±3 


125±6 











Table 1: Fits for different forms of \A\'^. We take row two as our result. 



Ignoring them in the fit does not affect the other parameters. Our final results for the Dalitz 
plot parameters are those shown in second row of the table. The corresponding correlation 
coefficients are shown in eq. (|8.2D . Fig |5| and Fig ^ show respectively a comparison between 
fit and data for the projections in X or y and the normalized residuals as function of bin 
number (left) and DP variables (right). 
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Figure 5: Comparison between data(points) and fit (histogram) for X,Y projections of tlie Daiitz 
plot distribution. 

normalized residuals , ,. , . , , , 




Figure 6: Left: Distribution of normalized residuals as function of bin number. Residuals fluctuate 
around zero, with 44 out of 154 exceeding 1, in absolute value. Right: Absolute value of normalized 
residuals distribution as function of X and Y . 

5. Systematic uncertainties 

We have estimated the systematics errors due to the following sources: 

Analysis cuts We have moved separately the following cuts: 6^^ in the range [15°, 21°] 
with a step of 3°, -P(x^) in the range [0.01, 0.15] with a step of 0.05, in the range 
[10, 25] MeV with a step of 5 MeV and J] in the range [780, 820] MeV with a step 
of 10 MeV. We find a negligible effect on the parameter estimates. 

Efficiency All reconstruction efficiencies have been checked with data, using control sam- 
ples. We find excellent agreement between data and MC for various kinematical 
distributions (see figj^. Concerning the photon detection efficiency we have checked 
that the error with which we estimate the ratio edata/^MC has a negligible impact 
on the estimate DP slope parameters. Only the EVCL procedure gives observable 
effects, as verified with the minimum bias sample. 

Resolution and binning Energy resolution for the photons is checked by comparing E^ 
distributions after the kinematic fit on data and MC. We find good agreement over 
the entire distributions. The drift chamber momentum resolution and absolute scale 
is checked run by run with the reconstructed Kg mass from Kg 7r"'"7r~ events. 
Binning size was changed up to a factor of two: 0.11 < AX, Ay < 0.2. 
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Background contamination The main source of backgrounds are: 4> ^ rj'j with r] — >■ 
vr^vr^vr'^ ,7r^ e^e^7 and (p — > lott^ with io 7r'''7r~7r'^ . Changing the cut on m^-y 
in a wide range, corresponding to a background change from 0.7% to 0.2% , we find 
smah changes for the parameter values. 

Stability with respect to data taking conditions We have divided our data sample 
in 9 periods of about 50 pb~^ each. We find that the results for each parameter are 
consistent with no change. 

Radiative corrections We have generated lO"^ r] tt^tt^tt^^ decays, according to ref. 
[0] . The bin by bin ratio of the DP density for r/ tt^tt^tt^ 7 decays to that for r/ — > 
vr+vr^vr'^ decays can be fitted with a constant with x^/dof = 154/153 corresponding 
to a CL of 46%. 

The results are shown in Table |2|. For each effect mentioned above the systematic error has 



Source 


Aa 


A6 


Ad 


A/ 


EVCL 


-0.017 


0.005 


-0.012 


0.01 


binning 


-0.008 +0.006 


-0.006 +0.006 


-0.007 +0.001 


-0.02 +0.02 


background 


-0.001 +0.006 


-0.008 +0.006 


-0.007 +0.007 


-0.01 


Total 


-0.019 +0.008 


± 0.010 


-0.016 +0.007 


± 0.02 



Table 2: Summary of the systematic errors on the Dalitz plot parameters. 

been estimated as the maximum parameter variation with respect to the reference value; 
the total systematic error is the sum in quadrature of the different contributions. 




100 200 300 400 500 -1 -0.8-0.6-0.4-0.2 -0 0.2 0.4 0.6 0.8 1 



Figure 7: Data vs Monte Carlo comparisons in log scale. Clockwise from top left: minimum pr 
and cosf? between pion tracks and i?-y for photons. 
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6. An alternative parametrization of the decay amplitude 



We have also fitted the Dahtz plot with a different parametrization which takes into account 
the final state tt-tt rescattering. Since strong interactions are expected to mix the two 
isospin /= 1 final states of the r/ — > Svr decay, it is possible to introduce a unique rescattering 
matrix R which mixes the corresponding 1=1 decay amplitudes [15| , for which we have: 



A 



A, 



(1) 



(1) / 

000 ^ R 



A 



(1) - 

+-0 



a: 



(1) 

000 



where: 



R = l + i 



a (5' 



(6.1) 



(6.2) 



(6.3) 



where oq = 0.18, Oq = —0.11, [3q = 0.06, /3q = —0.022 are obtained from ||16|] after proper 
rescaling from kaon to rj mass. The complete amplitudes, keeping the expansion in powers 
of X and y up to quadratic terms, are then given by: 



According to ref. [15|, the rescattering phases depend on the x and y variables^ as 

a = ao + 0{x^,y'^) (3 = (3q + 0{x,y) 
a' = a'o y + 0{x\y^) (3' = I3'^{y^ + x^/?,)/y + 0{x\ y^) 



(A+_o)fl = a{l + iao) - (6(1 + i/3o) + ia'^a) y + (c(l + iao) - d{l + i(3o) 
+i(3'ob) y^ + (c(l + iao) + d{l + i(3o) + x'^/3 



(6.4) 



and 



(^ooo)fl = 3 a( 1 + iao ) + [ 3 c(l + i ao) + 3i /3o 5 ]{x'^/3. + y'^) (6.5) 

We have fitted the Dalitz plot with the above parametrization and the fit results are given 
in Table ^. The systematic uncertainty on the parameters has been evaluated as described 



dof 




ax 1000 


6x1000 


cxlOOO 


dxlOOO 


150 


56% 


-7i.i2±o.o7l:[!:^^ 


13.71±0.04t°;°^ 


0.46±0.03l[!:J| 


-0.76±0.02t°:[J2 



Table 3: Results of the fit with a parametrization of the form eq. 6.4 



in section |5|. 

From the above results it is possible to extract the Dalitz plot slope a of the r] 
^0 ^0 ^0 (jgcay. From its definition: 



l^oooP oc 1 + 2az 



z = 9 mt/{4: ml Q^) x {x'^/3 + y"^) 



we get: 



a 



9 a (l + al) 



-0.038 ±0.003(stat);°;°J^(syst) (6.6) 



^We define x and y as: x — (si — S2)/m% and y = (ss — so)/m^ with Si = s,t,u for i — 1,2,3 and 
4 = {ml+ + ml- + m^o)/3. 
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in agreement with the PDG |^] average a = —0.031 it 0.004 and the recent KLOE prelim- 
inary result a = -0.027 ± 0.0041°;°°^ |l7|. 



7. Asymmetries 

While the polynomial fit of the Dalitz plot density gives valuable information on the matrix 
element, integrated asymmetries are very sensitive in assessing the possible presence of C 
violation in amplitudes of given A7. In particular left-right asymmetry - related to the c 
parameter in our fit - tests C violation with no specific AI constraint; quadrants asymmetry 
tests C violation for AI = 2 and sextants asymmetry (for a definition see ref. |]l8|) tests 
C violation for AI = 1. 

For this measurement care must be taken of possible slight differences between vr^ and 
7r~ reconstruction efficiencies. To this aim we estimate the MC efficiency separately for 
each region of the Dalitz plot, as the ratio between reconstructed and generated events in 
the region. This definition takes into account the resolution effects as well. From a sample 
of 5.7 X 10^ MC events we get: 

€L = (34.91±0.02)% CR = (35.05 ± 0.02)% 

ei3 = (35.01±0.02)% £24 = (34.95 ± 0.02)% (quad.) 

£135 = (35.00 ±0.02)% £246 = (34.96 ± 0.02)% (sext.) 
We have checked these values estimating the asymmetries on Monte Carlo: these turn out to 
be all compatible with zero. We then evaluate the asymmetries on data by subtracting the 
MC expected background and correcting the "raw" asymmetries with the above efficiencies. 
We obtain: 

^Lif = (9± 10)xlO"^ Aq = {-5 ±10) X 10-"^, = (8± 10)xlO"^ 

Systematic uncertainties on the asymmetries are obtained from studying: a) sensitivity to 
background, by varying cuts, b) event selection (EVCL) by use of the minimum bias sample 
and c) MC-data comparison using (p — > 7r"'"7r~7r^ events. In particular the tracking efficiency 
has been evaluated separately for the two charges, since in the MC a small but statistically 
significant difference in left and right efficiencies is evident. The above difference is due 
to a slightly different tracking efficiency vs pT for positive and negative pions because of 
nuclear interactions. 

Since we require both tracks to be reconstructed the absolute value of the efficiency 
is not important for the asymmetry, but rather its dependence upon the pion momentum. 
The good data-MC agreement has been already demonstrated for both charges on the 
signal. We here use the ip ~^ tt^tt^it^ control sample to check the agreement between data 
and MC for the tt"*" and vr" efficiencies as a function of momentum (see fig.^). 

The control sample agrees well with MC within errors, and the data-MC ratio is well 
fitted by a constant. 

In order to assess the possible systematic uncertainties connected with the tracking 
efficiencies we adopt a conservative approach: we estimat the maximum positive or negative 
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Figure 8: The data-MC ratio of tracking efficiency for 7r+ (top) and tt~ (bottom) vs pion pT- 



linear slopes compatible within one sigma with the fit of the distributions shown in fig]8|. 
Then we have assumed that the two charges behave with opposite slopes. This gives 
us two possibilities: vr"^ with positive slope and ir" with negative slope or vice-versa. 
We have then reweighted the events according to these two possibilities and used the 
maximum difference observed in the asymmetries as the corresponding systematic error. 
The systematics connected with the asymmetries are shown in Table 0. Therefore the 



Syst. Effect 


Left-Right 


Quadrant 


Sextant 


Background 
EVCL 
Efficiency 
Total 


(-0.2/ 0.1) X 10"^ 

(-0.5) X 10-3 
(-1.3/ + 0.9) X 10-3 
(-1.4/ + 0.9) X 10-3 


(-0.2/ + 0.2) X 10-3 

(-0.3) X 10-3 
(-0.3/ + 0.2) X 10-3 
(-0.5/ + 0.3) X 10-3 


(+0.3) X 10-3 
(+0.7) X 10-3 
(-1.3) X 10-3 
(-1.3/ + 0.8) X 10-3 



Table 4: Systematic errors on asymmetries, 
final results for the asymmetries are: 

^LR = (+0.09 + 0.10 1°:°^) X 10 
Aq = (-0.05 + 0.10 1°:°3) X 10 

= (+0.08 + 0.10 toil) X 10 
where the first (second) is the statistical (systematic) error. 

8. Conclusions 

The results including the statistical uncertainties coming from the fit and the estimate of 



-2 
-2 
-2 
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systematics are: 



a = -1.090 ± 0.005(stat)l°;°°°(syst) 
b = 0.124 ± 0.006(stat) ± O.OlO(syst) 
d = 0.057 ± 0.006(stat)l[]:[!?^(syst) 
/ = 0.14 ± O.Ol(stat) ± 0.02(syst) 



(8.1) 



Below we give the normalized correlation coefficients for the DP parameters. 



a 



b 



d 

-0.405 



/ 



a 1 

b 

d 



0.226 
1 



0.358 
1 



-0.795 
0.261 
0.113 



(8.2) 



/ 



1 



The following comments are in order: 

1. the fitted value for the quadratic slope in Y is almost one half of the simple current 
algebra prediction (6 = a^/4), thus calling for significant higher order corrections; 

2. the quadratic term in X is unambiguously different from zero; 

3. similarly for the large cubic term in Y; 

4. the fit results show correlations between the DP parameters. This should be properly 
taken into account for a correct error estimate when integrating the amplitude over 
phase space to get the decay width; 

5. fitting the r/ — > it'^it~it^ Dalitz plot with an alternative parametrisation we obtained 
a prediction for the r] 37r° slope which is consistent with the PDG average and 
the KLOE measurement; 

6. we do not observe any evidence for C violation in the rj — ^ 7r"'"7r~7r° decay since the 
c and e parameters of the Dalitz plot and the charge asymmetries are all perfectly 
consistent with zero. 
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